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ABSTRACT 

We study the response of the structure of an asymptotic giant branch (AGB) star 
to the accretion of a brown dwarf or planet in its interior. In particular, we examine 
the case in which the brown dwarf spirals-in, and the accreted matter is deposited 
at the base of the convective envelope and in the thin radiative shell surrounding 
the hydrogen burning shell. In our spherically symmetric simulations, we explore the 
effects of different accretion rates and we follow two scenarios in which the amounts 
of injected mass are equal to ~ 0.01 and ~ O.IMq. The calculations show that for 
high accretion rates (M acc = 10~ 4 Mq yr _1 ), the considerable release of accretion 
energy produces a substantial expansion of the star and gives rise to hot bottom 
burning at the base of the convective envelope. For somewhat lower accretion rates 
(M acc = 10~ 5 Mq yr _1 ), the accretion luminosity represents only a small fraction 
of the stellar luminosity, and as a result of the increase in mass (and concomitantly 
of the gravitational force), the star contracts. Our simulations also indicate that the 
triggering of thermal pulses is delayed (accelerated) if mass is injected at a slower 
(faster) rate. We analyze the effects of this accretion process on the surface chemical 
abundances and show that chemical modifications are mainly the result of deposition 
of fresh material rather than of active nucleosynthesis. Finally, we suggest that the 
accretion of brown dwarfs and planets can induce the ejection of shells around giant 
stars, increase their surface lithium abundance and lead to significant spin-up. The 
combination of these features is frequently observed among G and K giant stars. 

Key words: accretion - planetary systems - stars: AGB and post-AGB - stars: 
evolution - stars: low-mass, brown dwarfs - stars: abundances 



1 INTRODUCTION 

The discovery of several extra solar planets and brown dwarf 
candidates around solar type stars (e.g. Mayor & Queloz 
1995, Butler & Marcy 1996, Cochran & Hatzes 1997, Re- 
bolo et al. 1995, 1996, Basri et al. 1996) has advanced our 
understanding of the formation of planetary systems, but 
at the same time has also raised new astrophysical prob- 
lems. Notably, some of these planets are found to orbit very 
close to the central star (e.g. Mayor et al. 1997) and their 
orbits may be tidally unstable. For example, Rasio et al. 
(1996) showed that the orbital decay of 51 Peg's companion 
is inevitable (although on a long timescale), and the au- 
thors conclude that the planet will ultimately spiral-in into 
the star. Furthermore, the engulfing of planets is certainly 
a possibility when the star evolves to become a giant (e.g. 
Livio & Soker 1984). The present paper deals specifically 
with the reaction of a giant star "swallowing" a massive 
planet or a brown dwarf. This work was also motivated by 



the fact that the presence of a binary companion (even of a 
substellar mass) has often been invoked to explain the ax- 
isymmetrical morphology of elliptical and bipolar planetary 
nebulae (PNe) (e.g. Iben & Livio 1993; Soker 1992, 1996a; 
Livio 1997). Notably, if the substellar secondary is not too 
close to the primary, it will interact with the primary only 
on the upper AGB (e.g. Soker 1995, 1996b). 

In the course of its evolution, a main sequence star 
climbs along the red giant branch and eventually (at a later 
stage) can also reach the asymptotic giant branch (AGB). 
During the giant evolutionary phases, a planet may accrete 
from the wind emitted by the giant or be entirely engulfed 
by the giant's envelope. Depending on the initial orbital sep- 
aration, the planet or brown dwarf may be swallowed during 
the red giant phase, along the AGB, or survive. Once inside 
the envelope, due to viscous friction and also to tidal forces, 
a fraction of the orbital kinetic energy of the planet is con- 
verted into thermal energy and, as a consequence, the binary 
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separation decreases. During this process, orbital angular 
momentum is imparted to the envelope. The spiraling-in of 
a planet or a brown dwarf inside the envelope of an AGB 
star has been studied in several papers (e.g. Livio & Soker 
1984; Soker, Harpaz & Livio 1984). Using simplifying as- 
sumptions, these authors were able to follow the evolution 
of the brown dwarf's mass subject to accretion and mass 
loss (evaporation). Their model reveals the existence of a 
critical mass M cr i t below which the planet is evaporated or 
collides with the core. For the particular AGB star model 
used, they found M cr i t ~ 0.02Mq, but these results must be 
taken with caution since the physical processes involved in 
these calculations were treated only approximately. On the 
observational side, the orbital separations between close bi- 
nary stars seem to indicate that low-mass secondaries tend 
to have small orbital separations (e.g. Tables 2 and 3 in 
Iben & Livio 1993). A crude extrapolation of these data to 
low-mass stars (Soker 1996b) suggests that secondaries with 
masses of ^ 0.1 A/q will end up at orbital separations of 
IRq. As we will see in the next section, it is likely that at 
these depths inside giants stars, the secondary will be tidally 
disrupted or evaporated. Consequently, very low-mass sec- 
ondaries will probably be accreted onto the AGB star's core. 

In the present paper we follow the evolution of an AGB 
star that accretes matter from the destruction of a brown 
dwarf or a massive planet, using a spherically symmetric 
stellar evolution code. In the next section, we present the ini- 
tial model for the AGB star and the general input physics. 
In Section 3, we describe the numerical approach and in 
Section 4 we analyze the results of our computations. The 
effects of accretion on the nucleosynthesis and surface abun- 
dances are presented in Section H and a general discussion 
(Section |) and summary follow (Section M). 



2 INPUT PHYSICS AND VALUES OF 
PHYSICAL PARAMETERS 

To determine the exact location of the accretion process, we 
require a knowledge of the structure of the AGB star and the 
colliding brown dwarf. Where the brown dwarf will disrupt 
and dissipate can be estimated from the properties of the 
two objects involved in this scenario. In this section we thus 
present our initial model for the AGB star and look into 
the properties of the brown dwarf. With these data at hand, 
we will try to characterize the accretion process as well as 
different quantities associated with this phenomenon. 



2.1 The AGB star 

The initial AGB star model (taken from Forestini & Char- 
bonnel 1997) has an age of to = 4.74 10 s yr. Its initial ra- 
dius, mass, effective temperature and luminosity are equal 
to R = 198i? , M = 2.9M (initially M = 3.OM ), 
T e g = 3340K and L = 4395L , respectively. The star 
has already undergone 2 thermal pulses and is in the pro- 
cess of relaxing. The internal structure is composed of a 
very extended convective envelope which accounts for more 
than 99% of the radius. Below the massive convective zone 
(Menv — 2.36A/q), we find a thin radiative shell surrounding 
the hydrogen burning shell (HBS). In the HBS, H is mainly 
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burned through the CNO cycle. Below the HBS is the inter- 
shell where H is exhausted but where the temperature is too 
low to allow for helium burning. Below the intershell, helium 
is efficiently burned through 3a reactions in the so-called he- 
lium burning shell (HeBS). At the centre of the star, below 
the HeBS, the degenerate core is composed mainly of 16 
and 12 C. The core is highly degenerate and nearly isother- 
mal. In Table [j| we summarize the properties of the AGB 
star. Note that the maximum temperature occurs off-centre 
due to neutrino losses. 



2.2 The brown dwarf 

Brown dwarfs are compact, low mass objects (0.007AfQ 
Afbd 0.08A/q) which never undergo efficient hydrogen 
burning. Their evolution consists of a collapse phase, fol- 
lowed by deuterium burning (for the higher masses) and de- 
generate cooling (e.g. Stevenson 1991, Burrows et al. 1993). 
Zapolsky and Salpeter (1969) determined a mass- radius re- 
lation for this type of object, given approximately by the 
interpolation formula 
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Approximating the brown dwarf by a polytropic gaseous 
sphere with index n — 1.5 (characteristic of a non relativistic 
degenerate gas) , one derives for the central density, pressure 
and temperature the expressions 



Pc 



(3) 



(4) 



(5) 



where /ibd is the mean molecular weight, M aa \ and Rbd are 
the mass and radius of the brown dwarf, respectively. In Ta- 
ble [| we present several quantities associated with a brown 
dwarf and a very low-mass star of 0.015 and O.IMq, respec- 
tively. 

A straightforward comparison of the values in Tables 
[l] and | indicates that the brown dwarf's mean density is 
of the order of the density encountered in the HBS. It also 
reveals that brown dwarf radii are of the order of the radius 
of the hydrogen burning shell. 
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Table 2. Physical properties of accreting bodies 



M (M Q ) R(Rq) p bd (gem- 3 ) p c (gem" 3 ) P c (dynec m - 2 ) T c (K) 
0.015 0.095 24.6 147.6 2.4 10 16 1.9 10 6 

0.10 0.067 468.2 2806.2 4.3 10 18 1.8 10 7 



2.3 Times scales and orders of magnitude 

Given the properties of the AGB star and the brown dwarf, 
we can now attempt to determine the locus of the brown 
dwarf dissipation and the ensuing accretion process. 

First, we estimate the Virial temperature of the brown 
dwarf, since this represents the temperature above which the 
thermal kinetic energy of the star can overwhelm the grav- 
itational binding energy of the brown dwarf. This is given 
by 

Q Mbd Mbd m,H ,„ s 
Tv kR^— (6) 

Such temperatures, of the order of a few millions degrees, 
are encountered at the base of the convective envelope (see 
Table 0). Therefore, initially relatively massive brown dwarfs 
(M i; O.OIMq) will probably be able to reach the bottom of 
the convective envelope without being entirely evaporated. 
Note that we do not expect the trajectory of the spiraling- 
in brown dwarf to be strongly affected by the convective 
motions since its keplerian velocity is always supersonic (see 
e.g. Livio & Soker 1984). 

As the brown dwarf gets closer to the stellar core, tidal 
effects become important and they can induce strong dis- 
tortions to the brown dwarf's structure, and eventually they 
can even lead to a total disruption of the brown dwarf. The 
elongation stress at the brown dwarf centre can be approxi- 
mated by (e.g. Soker et al. 1987) 

„ f Rhd d (QM\ , Gp bi Rt d M _ 2 
S = J Phd dR[-W) ldl = ^ R? dyn6Cm ' 

(8) 

where pbd is the mean density of the brown dwarf, M and R 
the mass and radius of the star causing the disruption, and 
£ is a parameter of order unity. If we compare this expres- 
sion with the central pressure of the brown dwarf, then the 
condition of disruption, S > P c , reads 

( gp hd Rl d M \ V 3 / M \ V 3 / pbd \ l ' 3 

R < { P e ) ~ °' 27 Imq ) V 100 gem" 3 J 

/ R bi \ 2/ ' A ( Pc \-^ 3 R , , 
X IbTR^J ( lO^dynem-0 R& ' (9) 

From condition ^ , it is clear that a disruption of the brown 
dwarf occurs in regions close to the bottom of the convec- 
tive envelope. Thus, all the indications are that the brown 
dwarf will be dissipated (and its material accreted) near the 
bottom of the convective layer. 

In order to estimate the mean accretion rate onto the 
AGB star's core, we evaluate the orbital decay rate Tdecay in 
the accretion region. The spiraling-in of the brown dwarf is 
due to the action of drag forces that result from viscous fric- 
tion and tidal effects. This characteristic timescale, which is 
evaluated in the vicinity of the dissipation zone, provides a 



rough estimates of the transit timescale of the brown dwarf 
in this region and the rate at which matter can be deposited 
there. The accretion time scale is therefore given approxi- 
mately by (see e.g. Livio & Soker 1984 for a detailed discus- 
sion) 

- © 1 = 102 (a5r%) '(dy '" 

where p is the density in the envelope at radial distance r. 
Consequently, the expected accretion rate M acc is of the 
order of M acc ~ M bd /-r a cc ^ 6 10~ 5 M Q yr _1 for a O.OIMq 
brown dwarf. 

2.4 The accretion process 

The previous analysis indicates that a brown dwarf orbiting 
an AGB star can spiral-in to very deep layers inside the star. 
On the basis of our estimates of the Virial temperature and 
of tidal effects, we found that the brown dwarf is likely to be 
dissipated at the base of the convective envelope. The net 
outcome of the spiraling-in process is therefore the accretion 
of a brown dwarf mass at rates in the range of 10 to 

nn 4 Mqyt- 1 . 



3 NUMERICAL CALCULATIONS 

3.1 The stellar evolution code 

The stellar evolution code used in the present calcula- 
tions has been previously used for the computation of grids 
of standard pre-main sequence (PMS) evolutionary tracks 
(Forestini 1994, Siess et al. 1997a), for the detailed evolu- 
tion and nucleosynthesis of intermediate mass AGB stars 
(Forestini & Charbonnel 1997), and also for the effects of 
disc accretion onto PMS stars (Siess et al. 1997b, 1998). We 
refer the reader to these papers for a detailed description of 
the stellar evolution code. 

3.2 The models 

We computed two sets of models for two different accre- 
tion rates, M acc = 10~ 4 and 10~ 5 Mq yr" 1 (models A and 
B, respectively). The total accreted mass in cases A and B 
was O.O1245M and 0.01238A'/q , respectively. The chemical 
composition of the accreted matter was assumed to corre- 
spond to a solar metallicity (Z = 0.02) with relative abun- 
dances from Anders and Grevesse (1989). 

Previous works dealing with accretion in different evo- 
lutionary phases (e.g. Stahler 1988, Siess et al. 1997b) have 
pointed out the important role of deuterium in this process. 
Indeed, the very fast burning of 2 H through the 2 H(p,7) 3 He 
reaction releases a huge amount of energy which is able to 
modify the structure and evolution of the star significantly. 
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The presence or absence of deuterium in the accreted matter 
is thus an important issue which depends essentially on the 
mass of the brown dwarf. Brown dwarfs (or planets) more 
massive than ~ 12 Jupiter masses (~ O.OIMq, e.g. Bur- 
rows et al. 1995) have a high enough central temperature to 
burn deuterium in their interiors. Consequently, we do not 
expect deuterium to be accreted in the AGB star and its 
initial mass fraction has been set to be Xd = 10 -15 . We 
also assume that 7 Li is preserved in the brown dwarf and 
its mass fraction (abundance) is equal to Xu = 9.8968 10 -9 
[loge(Li)=3.31]. 



3.3 The deposition of accreted matter 



As we discussed in Section 2.4, a sufficiently massive brown 
dwarf is able to reach the bottom of the convective envelope. 
The deposition of some accreted matter inside the convective 
envelope, due to partial evaporation, will probably not in 
itself influence the structure since this region is devoid of 
energy sources. The main effects that need to be studied, 
therefore, are the consequences of mass deposition in the 
region close to the bottom of the convective zone, where the 
planet/brown dwarf is expected to disrupt. 

To achieve this goal, we have assumed a simplified mass 
deposition profile in which the mass AMj = M k +i — M k of 
the shell k is increased at every time step At by the amount 
fAMk, where / = M &cc At/(M tov - Mbot)- The mass is ac- 
creted in the region delineated by Mbot < M r < M top , where 
Mbot corresponds to the mass coordinate at the top of the 
HBS. However, if hot bottom burning is taking place (see be- 
low) , Mbot is equal to the bottom of the convective envelope 
(which is inside the HBS). With this prescription, we allow 
the accreted matter to go below the convective envelope, 
but, as we shall see later, a small fraction of the accreted 
mass (< 5%) is actually deposited in this region. The bound- 
ary Mtop is defined (arbitrarily) as the mass coordinate cor- 
responding to a radius equal to IRq (M top ~ 0.55Mq). 
This accounts for the fact that the planet is actually dissi- 
pated in the central region of the star. We also assume that 
the accreted matter is injected at the local envelope's tem- 
perature and density. This can be justified by the fact that 
the brown dwarf's kinetic energy is not too different from 
the specific thermal energy of the AGB envelope at the same 
point (e.g. Harpaz & Soker 1994). 

In the presence of mass accretion, at each evolutionary time- 
step, the total mass of the star is increased by the amount 
AM = M acc A£. In the shells where the mass is deposited 
(for Mbot < M r < Mtop), the time derivative of a variable 
A(M r ,t) is now calculated as 
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where u is the velocity of the shell and V = 
gravitational energy release, e gr av, is given by 



(11) 
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where ei n t is the specific internal energy. 

The chemical composition entering the computation at 
the new time step is modified due to the arrival of fresh, 
not nuclearly processed material from the accreted matter. 



The mass fraction X(t + At)i,k of element i at shell k after 
mixing, is given by (Siess et al. 1997b) 



X(t + At) i>k = 



X(t)j, k + fX? 

1 + / 



(13) 



where X(t)i t k is the stellar mass fraction of element i be- 
fore mixing and X" is the mass fraction this element in the 
accreted matter. Note that X" is constant and does not de- 
pend on the mass coordinate. 

Finally, in order to speed-up the convergence process, 
we turned off the wind mass loss during the accretion phase 
and we used grey atmosphere models. 



4 EVOLUTION OF THE STRUCTURE 

The evolution of the AGB star under consideration can be 
divided into 2 major phases. The first one deals with the 
onset of accretion and the main accretion phase while the 
second marks the end of the accretion process and the re- 
laxation of the star. 



4.1 The onset of accretion and the main accretion 
phase 

The deposition of material inside the star releases a signif- 
icant amount of gravitational energy. The associated accre- 
tion luminosity is given by 



GMM a 
R 



M a , 



lO- 4 M yr- 1 
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(14) 



Note that an additional source of energy can arise from the 
possible replenishment of the nuclearly active region with 
accreted material. The modification of the chemical compo- 
sition of these regions can increase the nuclear energy pro- 
duction. 

For high accretion rates (case A), L acc represents a large 
fraction of the stellar luminosity and in the deep region 
where the energy is deposited, the star expands. Locally, 
the density, pressure and temperature decrease. The HBS, 
which advances in mass, rapidly reaches the "cooler" accre- 
tion region^] and then, the nuclear energy production due 
to H burning decreases (see Fig. [j]). On the other hand, 
the luminosity of the HeBS increases, since the addition 
of mass results in an increased gravity and compressional 
heating. In the expanding accretion region, the opacity k 
increases since in this regime it is a decreasing function 
of the temperature. The increase in both k and L (due to 
Lacc) contributes to a steepening of the radiative gradient 
(V ra d oc kL) which now becomes greater than the adia- 
batic gradient, thus convection sets in. Consequently, the 
convective envelope penetrates rapidly towards deeper re- 
gions (see Fig. ^) and reaches the HBS, leading to the so- 
called hot bottom burning (HBB). With the onset of HBB, 



For the following discussion, we identify the accretion region as 
the region located above the HBS and up to a mass coordinate 
M r ~ 0.55Mq where the mass is deposited. Note also that the 
top of the accretion region is located in the convective envelope. 
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Figure 1. The evolution of the radius (R), nuclear luminosity 
(L nuc ) and nuclear luminosities associated with H (Lh) and He 
(Ljje) burning, in the case of high accretion rate (case A). Open 
and filled dots mark the beginning of the hot bottom burning 
and the end of the accretion process, respectively. In case A, 
0.01245Mq nas been accreted in the star. 



the luminosity of the HBS (Lh) rises sharply, as fresh ma- 
terial coming from the convective region is engulfed in this 
shell. The global contribution to the total nuclear energy 
production of the 13 C(p,7) 14 N reaction is increased during 
this phase, and it accounts for 27% (compared to 22% in 
a standard evolution). We also note that the 14 N(p,7) 15 
and 15 N(p,a) 16 reactions contribute less to the nuclear 
energy production than in a standard scheme. These dif- 
ferent features indicate that the CNO cycle is reaching a 
new equilibrium. The "flash" induced by the penetration 
of the convective envelope into the HBS increases Lh by 4 
orders of magnitudes. The large release of nuclear energy 
is used to expand the envelope and shortly after the igni- 
tion of the HBB, Lh decreases rapidly (Fig. |l|). Thereafter, 
La.cc provides most of the luminosity, the HBS almost extin- 
guishes as it moves into the colder accretion region, whereas 
Lhc rises since the mass increases. Figure ^ depicts the pro- 
files of several physical quantities during the early accretion 
phase. The peaks in the luminosity profiles indicate that the 
star is not in thermal equilibrium, but rather energy accu- 
mulates in its interior. One also notices that this feature is 
smoothed out as time advances, as a result of the star's ex- 
pansion. In our initial model (solid line), the gravitational 
energy production £ gra v represents a small fraction of the 
stellar energy production, 80% of the luminosity is due to 
nuclear energy. When the accretion process starts (dotted 
line), £ g rav increases rapidly below the accretion region as 
a result of mass deposition. In the accretion region where 
£grav < 0, the expansion induces modifications to the tem- 
perature profile which is now steeper. The nuclear energy 
production is consequently more localized and concentrated 
to thinner regions. When HBB starts (short-dashed line), 



Figure 2. The evolution of the boundaries of the HeBS, HBS 
and convective envelope. The dotted lines indicate the locus of 
maximum energy generation while the dashed lines delineate the 
burning shells (e n uc > 5 ergg - ^^s - 1 ). Note that the bottom of 
the HeBS has not been marked in the graphs, its location is at 
about M r ~ 0.47Mq. The solid line indicates the lower edge 
of the convective envelope. The evolution is presented for the 
standard case (left panel) and for cases A and B (middle and right 
panels, respectively). The important feature is the appearance 
of hot bottom burning in case A. The top panel represents an 
enlarged view of the HBS in case A. 

both Egrav and e nU c reach extremum values. The luminos- 
ity becomes negative, as most of the energy production is 
used to expand the star. Finally, when the HBB is settled 
(long-dashed line) , the star heats up in the region of nuclear 
energy production (e gra v > 0), whereas expansion proceeds 
above this zone. During the accretion phase the HBS nar- 
rows as the temperature profile becomes steeper (Fig. 
This also accounts for the global decrease of L nuc over this 
period. Finally, since mass is being added, all the structural 
features (location of the maximum energy production, burn- 
ing shells, step in temperature) move towards higher mass 
coordinates. 

In case B, the release of accretion energy is not sufficient 
to expand the star and trigger convective energy transport 
down to the HBS. Consequently, a thin radiative layer is 
maintained below the convective envelope, like in the stan- 
dard case (Fig. ^). This thin layer prevents the energy pro- 
duction of the HBS from heating the envelope directly and, 
as L acc is small and the gravitational pull increases, the star 
contracts. At the beginning of the accretion process, the re- 
lease of gravitational energy lowers the temperature in the 
accretion region (above the HBS). As the core grows due to 
mass addition, the HBS moves outwards and rapidly pene- 
trates the less dense and cooler accretion region. As a result, 
the nuclear luminosity due to hydrogen burning (Lh) under- 
goes a sudden drop (Fig. ^J). The star compensates for this 
decrease in energy production by a small contraction. There- 
after, the star tries to find an equilibrium configuration in 
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Figure 3. The evolution in the case of high accretion rates of the 
profiles of gravitational (nuclear) energy production rates £g ra v 
(£mic)> temperature (T) and luminosity L r in the accretion re- 
gion. The top panel shows the profile of L r over the entire star. 
The solid curve refers to our initial model, without accretion. 
The dotted, short- and long-dashed lines correspond to models 
for which accretion is taking place for 15, 27 and 67 years, re- 
spectively. 



which L nU c and L acc account for the total luminosity. This 
phase of structural readjustment leads to a small expansion. 
Thermal equilibrium is finally achieved after ~ 130 yr, which 
corresponds to the stellar Kelvin Helmholtz timescale 
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yr 



(15) 



Subsequently, contraction resumes and the nuclear luminos- 
ity stays at an almost constant value, somewhat lower than 
in the standard case (without accretion), due to the pres- 
ence of L a cc. However, the increase in mass speeds-up the 
contraction rate in the central region and the luminosity of 
the HeBS rises. 

The main difference between case B and the standard 
evolution is the persistence of stellar contraction, not only 
in the core but throughout the entire structure during most 
of the accretion phase. For low accretion rates, it is mainly 
the mechanical effects of mass deposition and gravitational 
attraction that prevail. For a comparison with the high ac- 
cretion rate case, we present in Fig. ^ the profiles of several 
physical variables during this period. As can be seen, e gra v 
remains largely positive in the accretion region : contrac- 
tion heats the central parts of the star. The profiles of T, 
Enuc and L are very similar during the entire phase. We 
do not obtain a significant increase in the steepness of the 
temperature gradient nor a peak in the luminosity profile. 
The different quantities are simply shifted to higher mass 
coordinates very smoothly. In this case, accretion does not 



Figure 4. The evolution of the radius (R), nuclear luminosity 
(Lnuc) and nuclear energy production due to H (Ljj) and He 
(Ljje) burning. The solid line refers to case B, the dotted line to 
the standard model. Filled circles mark the end of the accretion 
phase. 




L r (L G ) - 3000 



Figure 5. The evolution (in the case of low accretion rates) of 
the profiles of luminosity (L r ), gravitational (nuclear) energy pro- 
duction rates e grav (e n uc) and temperature (T) in the accretion 
region. The solid curves refer to our initial model, without accre- 
tion. The dotted and short-dashed lines correspond to models for 
which accretion is taking place for 327 and 535 years, respectively. 
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Figure 6. Relaxation of the star. The variables are the same as 
in Fig. ^. The solid, dotted and dashed lines refer to cases A, B 
and to the standard evolution models, respectively. 



lead to significant structural changes, the energetic effects 
of accretion (L acc ) being too weak. 



4.2 The relaxation of the star and the subsequent 
evolution 

For high accretion rates (Fig. [l]), the sudden suppression of 
accretion energy deposition initiates a rapid contraction of 
the star, as the envelope can no longer be supported. The 
release of gravitational energy increases the temperature in 
the regions of energy production and induces an enhance- 
ment of the nuclear luminosities. The star globally heats 
up, the opacity drops, and the convective envelope recedes 
towards the surface (Fig. The cessation of the accretion 
process thus marks the end of HBB and this leads to the 
formation a thin radiative layer between the convective zone 
and the HBS. The contraction proceeds as long as the en- 
ergy losses from the surface are not equilibrated by the gain 
from gravitational and nuclear energy sources. The return to 
equilibrium is finally achieved after ~ 675 yr, when the stel- 
lar radius and luminosity have decreased by a factor of ~ 2 
(Fig. ^|). At this point, the nuclear luminosity provides more 
than 90% of the stellar luminosity and it is sufficient to sup- 
port the envelope. As most of the nuclear energy is produced 
by the HBS, the energetic influence of the HeBS is shielded 
and its luminosity (Lh c ) decreases (as in the standard case). 
Finally, as the core mass M core grows, the temperature in- 
side the HeBS increases and the nuclear reactions (which are 
very temperature sensitive) runaway and trigger a thermal 
instability. Note that in case A, the thermal pulse is delayed 
with respect to the standard model. This is due to the large 
expansion that occurs during the accretion phase, which in 
turn slows down the "natural" temperature increase of the 
HeBS substantially. 

In case B, when the accretion process stops, the star 



undergoes a short phase of rapid contraction with the dis- 
appearance of L acc . The temperature in the HBS increases 
and the nuclear luminosity rises (Fig. ^) . For a short period 
of time the star expands. Finally, when an equilibrium con- 
figuration is reached, which accounts for the missing L acc , 
the star resumes contraction and LHe keeps rising. The ten- 
dency of increasing Lhc is due to the fact that for low accre- 
tion rates the impact of accretion is more a consequence 
of the mechanical effects resulting from mass deposition 
(compressional heating), rather than energetic effects due 
to L acc or HBB. Note in addition that the accretion rate 
used in this study is larger that the "normal" core growth 
rate (~ 10 _7 MQyr _1 , for a 3A/q star). Therefore, accre- 
tion forces the star to evolve more rapidly than in a standard 
evolution and the thermal pulse occurs earlier. However, the 
thermal instability is weaker in the star accreting the brown 
dwarf than in the standard model: The maximum of Lh c 
is more than an order of magnitude lower than in the con- 
ventional computation. The reason for this difference is the 
fact that due to the mass deposition, the HeBS moves to- 
wards regions of a lower density faster than in the standard 
scheme, and the nuclear energy generation is thus reduced. 

4.3 Accretion of more massive bodies 

If a high accretion rate (M acc = 10 -4 Mq yr -1 ) is sustained 
for a longer period of time, the release of gravitational en- 
ergy maintains the envelope expansion. After ~ 240 yr, the 
Kelvin Helmholtz timescale has decreased to ~ 60 yr and 
accretion energy can now be efficiently radiated away. At 
that time, the nuclear energy production is at a minimum 
and the HBS is almost extinguished (Fig. [?]). The luminos- 
ity profile is now monotonicly increasing and contraction 
resumes. When the accretion process ends (for an accreted 
mass M acc = 0.0594Mq in this simulation), the contrac- 
tion rate accelerates as the envelope is no longer supported 
by the deposition of accretion energy. The temperature in- 
creases and LHe rises first, followed by Lh. When the nu- 
clear sources provide Jj 90% of the stellar luminosity the 
contraction stops. The radius is then minimal and its value 
is smaller than in case A, because in this simulation, the 
stellar mass and thus the gravitational pull are larger. Fi- 
nally, as the temperature increases the thermal instability 
sets in, but somewhat later than in case A. 

For lower accretion rates (M acc = 10~ 5 Mq yr -1 ), the 
persistent addition of mass, up to M acc = 0.1082Mq, in- 
creases the contraction rate. The star heats up faster and 
the thermal instability develops even earlier than in case B 
(Fi g. 8). For the reasons explained previously (see Section 
4.2), the value of Lh c during the pulse is lower than in the 
standard case. We also note that accretion does not perturb 
the thermal pulse, which indicates that in this case this pro- 
cess represents a small perturbation to the stellar structure. 
When the instability is over, the star resumes contraction. 
The nuclear luminosity (Z/ nuc ) is lower than in the standard 
model because a fraction of the energy is provided by the 
accretion process. Finally, when accretion ends, the star re- 
laxes and returns to a standard evolution (with a larger mass 
however) . 

To summarize these simulations, for high accretion rates 
the accretion luminosity dominates the stellar energetics, re- 
ducing the nuclear energy production of the HBS. As a re- 
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Figure 7. The influence of mass deposition for high accretion 
rates (10 — 4 Mq yr — 1 ). The physical quantities are the same as 
m Fig. M. The solid curves correspond to a long-lasting accretion 
phase (Mace = 0.0594MQ), the dotted line to case A (M acc = 
O.O1245M0) and the dashed line is the standard model. 



Figure 8. Same as Fig. Q but for low accretion rates 
(10~ 5 Mq yr _1 ). The solid curves correspond to M acc = 
0.1082Mq, the dotted lines to case B (M acc = 0.01238M Q ) and 
the dashed lines represent the standard model. 



suit, the (nuclear) evolution of the star is slowed down. For 
lower accretion rates (but still high compared to the normal 
core growth rate), accretion shifts the internal structure to 
that corresponding to higher masses at a faster rate than 
in a standard evolution, and consequently the evolution is 
accelerated. However, these results do not imply that for 
some intermediate value of M acc the evolution of the star re- 
mains unchanged. Indeed, the stellar mass increases and the 
core's growth rate remains always higher than in a stan- 
dard evolution. In the limiting case where the release of 
accretion energy exactly counterbalances the gravitational 
attraction, instantaneously the evolution of the star would 
remain the same. However this situation is unstable since 
M core increases. Quite rapidly the gravitational pull takes 
over and the contraction resumes. Note that for large ac- 
cretion rates when the accretion process is maintained for a 
longer period of time, contraction finally resumes once ther- 
mal equilibrium is reached. 



5 SURFACE ABUNDANCES AND 
NUCLEOSYNTHESIS 

Chemical modifications resulting from mass accretion are 
expected if either a substantial amount of accreted matter 
is deposited in the envelope or/and if nuclear burning is 
operating in the convective envelope. To separate these two 
effects, we define new variables Ri and R2 for each chemical 
element i by 



R 2 ,i = 



Xi(t) 



Ri 



Xj(t) 

*?(*) 



XI{t) + mMoX° 



(17) 



1 + m a 



(16) 



where Xi(t) is the mass fraction of element i in an evolu- 
tionary scheme with accretion and X? (t) is the mass fraction 
of this element in standard computations. m acc represents 
the ratio of the total mass deposited in the envelope during 
the accretion process (Af acCiCnv ) to the mass of this region 
(Afenv). The variable Ri describes the effect of accretion on 
the surface abundances, including mass deposition and nu- 
clear burning. The variable R2 represents the nuclear activ- 
ity: it gives the ratio of the mass fraction computed in our 
models with accretion to the mass fraction expected from 
calculations with accretion if no nuclear burning took place. 
Therefore, if the element is nuclearly produced then i?2 > 1 
and if it is destroyed R2 < 1. 

The deposition of fresh, nuclearly unaltered material 
in the convective envelope entails modifications of the sur- 
face chemical abundances. This effect accounts for the large 
values of Ri attached to light elements such as 6 Li, 7 Li, 
9 Be, 10 B or n B (Fig. ||). Indeed, due to their low burning 
temperature, these elements were previously burned dur- 
ing the pre-main sequence phase or during the first dredge 
up. Depending on the mass added to the envelope, the 
abundances of these elements can be increased by a fac- 
tor of ~ 2 — 8 (0.3 — 0.9 dex). The more mass is accreted, 
the higher the ratio Ri is. However, the detection of these 
chemical changes may be difficult if the accreted mass is 
small. For example, the 7 Li mass abundance in cases A 
and B (M acCj0nv ~ O.OIMq) reaches values of log e(Li) = 1.54 
and log e(Li) = 1.51 respectively, which are very close to the 
log e(Li)=1.40 standard value. It is interesting to note that 
the accretion of planetary material in the envelope of red 
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Figure 9. The effects of accretion on the surface abundances of 
several chemical elements. The lower panel shows the ratio of the 
mass fraction of an element resulting from our computations with 
accretion to the standard value. The solid and short-dashed lines 
correspond to cases A and B, respectively. The dotted and long- 
dashed lines correspond to longer accretion phase, with accretion 
rates (accreted masses) that are equal to 10 — 4 (O.O594M0) and 
10 — 5 Mq yr~ 1 (0.1082Mq), respectively. The upper panel repre- 
sents the nuclear activity. If R2 > 1 (R2 < 1) the element is 
nuclear ly produced (destroyed, see text). 

giant stars was originally suggested by Alexander (1967) 
to explain the high Li abundance observed in some red gi- 
ants. Brown et al. (1989) showed that a sufficient mass of 
terrestrial-like planets (which have lost their initial H and 
He content) may replenish the envelope of giants sufficiently 
to the reproduce high observed Li abundances. Note, fi- 
nally, that the surface abundances of heavier elements, with 
A > 13, are close to their initial values and are thus not 
modified by the accretion process. 

The upper panel of Fig. |^ indicates which elements un- 
dergo efficient nuclear reactions. The light elements are the 
principal target of nuclear burning. Indeed the "low" tem- 
perature at the base of the convective region, between 3 10 6 
and 10 7 K (depending on the presence or the absence of 
HBB), allows the destruction of these species. The nuclear 
burning acts preferentially on 6 Li, 7 Li and 9 Be. We can also 
detect the effects of CNO reactions in the envelope through 
the production of 13 C, 14 N and 17 and the depletion of 12 C. 
However, the influence of H burning through the CNO cycle 
is very small and it only appears if accretion is maintained 
for a sufficiently long period of time (dotted line, dashed line 
in Fig. g). 

To summarize, the temperature at the base of the con- 
vective envelope (Ss 10 7 K) is too low to give rise to a rich 
nucleosynthesis. The main chemical changes thus result from 
the deposition of fresh material in the convective envelope. 



6 DISCUSSION 

Our computations involved a relatively young AGB star (at 
the beginning of the thermally pulsing phase) and one may 
wonder what would be the effects of brown dwarf accretion 
inside more evolved stars. In this respect we note, that as 
the evolution progresses, the core grows but the temperature 
at the base of the convective zone remains almost constant, 
close to 4 10 6 K. Thus, the location of the accretion process 
may remain unchanged. However, in older AGB stars the 
stellar luminosity is substantially higher and this will most 
certainly result in a different stellar response to high ac- 
cretion rates. Indeed, as we have seen in the present work, 
the main reason for the stellar expansion is the release of 
a large amount of accretion energy (L acc > L»). We thus 
expect that if mass is injected in more luminous AGB stars, 
then the stellar radius and luminosity will decrease during 
the process, as in case B (see also Harpaz & Soker 1994). 

Our calculations revealed that for high mass injection 
rates (M acc = 10 — AfQyr -1 ) the radius and luminosity of 
the star increase by a factor of ~ 2. This leads to a substan- 
tial increase in the mass loss rate, since according to the 
Reimers formula (1975) Mi oss cc LR/M. The ejected mass 
can be estimated as 




where £ aC c is the duration of the accretion process. We can 
therefore expect the formation of shells, corresponding to 
individual accretion events. The mass of these shells is of 
the order of ~ 1(T 3 - 1O" 4 M . 

A phenomenon not directly addressed by our spherically 
symmetric calculations is that of the deposition of angular 
momentum into the giant's envelope, by the accreted brown 
dwarf. We can however estimate the effects of such a depo- 
sition as follows. If we assume that the secondary deposits 
all of its orbital angular momentum in the envelope, and we 
further assume that the envelope rotates rigidly (due to the 
fact that it is convective, and thus viscous) then the ratio 
of the envelope angular velocity ui to the critical velocity 
(surface keplerian angular velocity) luk is given by 

^■""•(iiisfeXsf ■ W 

where a is the initial orbital radius of the brown dwarf. Here 
we calculated the moment of inertia of the envelope for our 
model, and it is given by J cnv = k 2 M ellv R 2 , where k 2 ~ 
0.12. As we can see, Eq. y_9j indicates that for typical brown 
dwarf masses, the AGB star can be spun-up to a significant 
fraction of the critical velocity. Notably, the accretion of 
angular momentum is required to explain the abnormally 
high rotational velocities of FK Comae stars (e.g. Rucinski 
1990, Welty & Ramsey 1994) and of peculiar AGB stars such 
as V Hydrae (Barnbaum et al. 1995). Finally, this can also 
have important consequences for the shaping of planetary 
nebulae (see e.g. Soker 1996a). In particular, it has been 
suggested that the spin-up of the envelope is responsible 
for the high fraction of non-spherical planetary nebulae (see 
Livio 1997 for a review). 

Concerning peculiar abundances, the discovery of 
lithium rich G and K giant stars (e.g. Wallerstein & Sne- 
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den 1982, Brown et al. 1989, de la Reza & da Silva 1995) 
posed the problem of the origin of this element in these stars. 
Indeed, standard stellar evolution predicts that Li should be 
depleted during the pre-main sequence phase or during the 
first dredge up. Several mechanisms have been proposed to 
explain the high Li abundance, these include : (1) Li pro- 
duction (Cameron & Fowler 1971), (2) the retainment of the 
original lithium, (3) the external enrichment by novae (e.g. 
Gratton & D'Antona 1989), (4) the engulfing of a planet 
(e.g. Brown et al. 1989). Recently, de la Reza et al. (1996, 
1997) pointed out that almost all the Li K giants are optical 
counterparts of IRAS sources, the far IR excesses being at- 
tributed to the presence of a dusty circumstellar shell (CS). 
Based on this apparent correlation between the IR and Li 
excesses, these authors proposed a scenario in which all the 
K giants with masses between 1.0 and 2.5 Mq become Li 
rich after the first dredge-up. They assume that the mecha- 
nism responsible for the Li enrichment is accompanied by a 
sudden mass loss which leads to the formation of a CS. The 
Li which has not been ejected remains in the stellar surface 
and is depleted later when the mass loss stops. In this sce- 
nario, the CS expands, cools down, disappears and the Li 
surface abundance reaches again standard values. This evo- 
lutionary path, or "Li path", can reasonably explain the Li 
abundance distribution of K giant stars in a colour-colour 
diagram. Its duration is of the order of ~ 10 yr and it re- 
quires that Li suddenly appear at the surface of the star 
and then be depleted on a timescale of 10 3 — 10 5 yr. How- 
ever, de la Reza et al. do not provide convincing mechanisms 
for Li enrichment and mass loss. They invoke an extra mix- 
ing process for Li production, the "cool bottom processing" 
(Wasserburg, Boothroyd & Sackmann 1995), but they re- 
main vague concerning the generation of mass loss. We pro- 
pose here that the accretion of a planet or a brown dwarf 
by a giant star could effectively explain all of these features. 
Indeed, our calculations show that the accretion process (1) 
can increase the mass loss rate and lead to the ejection of a 
shell (Eq. [jj), and (2) can enrich the envelope with lithium 
(Fig. H), in agreement with observations. 

Observations also show that lithium rich giants have 
generally normal rotational velocities [de Medeiros et al. 
1996; although some rapid rotators exist, e.g. IE 1751+7046 
(Ambruster et al. 1997)]. This implies that if the accretion 
of a planet/brown dwarf is responsible for the Li enhance- 
ment, some fraction of the orbital angular momentum which 
is deposited into the envelope has to be removed (or that in 
most cases of engulfed planets a/R < 1, see Eq. |l9|). This can 
be achieved by efficient magnetic braking, since the giants 
develop a deep convective envelope (e.g. Leonard & Livio 
1995). We should note that the paucity of giants showing 
both rapid rotation and high Li abundance probably indi- 
cates that the timescale for rotational braking is short com- 
pared to the giant evolutionary phase (a few 10 7 yr). 

The existence of very rich Li stars (e.g. Gregorio-Hetem 
et al. 1992, Torres et al. 1995, de la Reza & da Silva 1995, 
Fekel et al. 1996) poses a somewhat more difficult challenge. 
One could imagine, however, that the collision of the planet 
(or a brown dwarf) with the star's core may trigger the "cool 
bottom processing" (Wasserburg, Boothroyd & Sackmann 
1995). In this scenario, 7 Be-rich material can be dredged-up 
from the deep interior to the convective envelope where it 
decays through the 7 Be(e~, ^) 7 Li reaction. 



A comment should be added about statistics. From a 
sample of 644 stars, Brown et al. (1989) found that ~ 4% 
of G and K giants have Li abundances above log e(Li) ~ 1.3 
and that an additional 4% have abundances between 1.2 and 
1.3. Considering the fact that many of the Li rich stars show 
an IR excess compatible with the presence of a surrounding 
shell, it is tempting to interpret these observations as a sig- 
nature of an accretion event. This could indicate that 4 — 8% 
of all stars have low-mass (possibly substellar) companions. 
These statistics are fully compatible with those inferred from 
PNe nuclei (e.g. Yungelson, Tutukov & Livio 1993). 



7 SUMMARY 

We have presented computations in which a brown dwarf (or 
a planet) is accreted in the core of an AGB star. This situa- 
tion results when the brown dwarf is engulfed in the stellar 
envelope as the latter is expanding due to stellar evolution. 
Viscous friction as well as tidal effects generate a drag force 
that makes the brown dwarf spiral in. We have shown that 
the region where the brown dwarf is expected to disrupt is 
close to the base of the convective envelope. 

We investigated the response of the star to the ac- 
cretion of a small body for two accretion rates, 10 -4 and 
10~ 5 Mq yr -1 . Our simulations show that for the higher ac- 
cretion rate, the main energetic effects come from the release 
of potential energy. The accretion luminosity is of the order 
of the stellar luminosity and this additional energy source 
makes the star expand. Another interesting feature of these 
computations (involving a high accretion rate) is the ap- 
pearance of hot bottom burning. However, the relatively low 
temperature at the base of the convective envelope (~ 10 7 K) 
and the short duration of the accretion phase do not allow 
for a rich nucleosynthesis in the convective zone. We find 
that accretion delays the development of the thermal insta- 
bility. The nuclear energy production due to the H and He 
burning shells is smaller than in a standard case because a 
part of the energy is provided by the deposition of the accre- 
tion energy. Hydrogen and helium are consequently burned 
more slowly and the evolution is slowed down. 

For the lower accretion rate the evolution is different. In 
this case the accretion luminosity represents a small fraction 
of the stellar luminosity and therefore expansion does not 
ensue. Conversely, the addition of mass increases the grav- 
itational pull and the star contracts even more efficiently 
than in a standard scheme. The temperature of the nuclear 
burning shells increases more rapidly and the thermal in- 
stability occurs earlier than in conventional evolution. The 
global stellar evolution is accelerated. 

Except for variations in the luminosity and radius, both 
of which can be increased by a factor of ~ 2, the effects of 
mass accretion in an evolved star may be traced by changes 
in the surface chemical abundances. We showed that such 
changes are mainly the result of material deposition rather 
than nucleosynthesis. The chemical changes affect mainly 
light elements such as 7 Li and 9 Be whose mass fraction can 
be increased by a factor of 2-8, depending on the accretion 
rate and duration of the accretion process. 

Generally, the observational signatures of the accretion 
of a brown dwarf or a planet by giant stars can be (depending 
on the physical parameters involved): (i) the formation of 
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a circumstellar shell, (it) the enhancement of the surface 
7 Li abundance, (Hi) rapid rotation. We propose that the 
IR excesses and high Li abundances observed in 4 — 8% of 
the G and K giants originate from the accretion of massive 
planets/brown dwarfs/very low-mass stars. 

Finally, we should note that we also made several at- 
tempts to accrete matter deeper inside the star (e.g. below 
the HBS). However, such a process could not be followed re- 
liably with our code. The deposition of protons in the very 
central regions increases the nuclear energy production to 
such rates, that the envelope cannot absorb this energy ex- 
cess. 
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